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We report a neutron scattering study of the long-wavelength dynamic spin correlations in the model 
two-dimensional S = 1/2 square lattice Heisenberg antiferromagnets Sr2Cu02Cl2 and Sr^ChiaCuCla. 
The characteristic energy scale, Uo(T/J), is determined by measuring the quasielastic peak width in 
the paramagnetic phase over a wide range of temperature (0.2 < T/J < 0.7). The obtained values 
for uiq[T / J) agree quantitatively between the two compounds and also with values deduced from 
quantum Monte Carlo simulations. The combined data show scaling behavior, u ~ , over the 
entire temperature range with z = 1.0(1), in agreement with dynamic scaling theory. 

PACS numbers: 75.10.Jm, 75.40.Gb, 74.72.-h 



Since the discovery of high temperature superconduc- 
tivity in 1986 a great number of studies have been de- 
voted to understanding this fascinating yet difficult prob- 
lem. Although the superconducting mechanism itself in 
these lamellar copper oxides is still elusive, tremendous 
progress in condensed matter physics has resulted as a 
byproduct. One such example is our understanding of 
the two-dimensional quantum Heisenberg antiferromag- 
net (2DQHA) on a square lattice, which had been a long- 
standing problem in theoretical physics even before the 
discovery of high temperature superconductivity. Since 
the parent compounds of the copper oxide superconduc- 
tors, such as La2Cu04, are very good representations of 
the 5 = 1/2 square-lattice Heisenberg antiferromagnet, 
the study of the 2DQHA has received renewed interest 
from both theorists and experimentalists fjjj. The syner- 
gistic efforts of theoretical, numerical, and experimental 
investigations have yielded a comprehensive picture of 
the static magnetic properties of the 2DQHA over the 
past decade 

However, less progress has been made in studies of 
the dynamic properties of the 2DQHA. In particular, 
there are only a small number of analytic theoretical 
models and numerical simulations on this problem, while 
systematic experimental investigation is lacking for the 
5 = 1/2 2DQHA. One of the most fundamental questions 
we can ask in the study of dynamic critical behavior of 
the 2DQHA is whether or not dynamic scaling is obeyed. 
The dynamic scaling hypothesis proposed by Ferrell et 
al. Q| and developed for magnetic systems by Halperin 
and Hohenberg |3|, is that the dynamic structure factor 
is completely determined by the static properties with 
the scaling frequency given by u>q ~ £,~ z - The theo- 
retical value for the dynamic critical exponent z = d/2 
of the Heisenberg antiferromagnet has been known for 
many years 0. Using this dynamic scaling hypothesis, 



Chakravarty, Halperin, and Nelson (CHN) have argued 
that in the 5=1/2 2DQHA the characteristic energy 
scale is given by uj cx £- 1 c {T/2irp s ) 1 / 2 , where Co is the 
spin-wave velocity and p° s is the spin stiffness constant, 
both at zero temperature We set K = Hb = 1, so 
that temperature T, and frequency u> have the units of 
energy. 

The characteristic energy scale, ojo, for the 2DQHA is 
equivalent to the relaxation rate of the order parameter, 
which can be determined by measuring the damping of 
the quasielastic peak at the antiferromagnetic wave vec- 
tor (7T, 7r) via inelastic neutron scattering. One of the 
main reasons for the scarcity of neutron experiments is 
the large energy scale of the magnetism in the copper ox- 
ide materials. In magnetic systems, the primary energy 
scale is set by the nearest neighbor exchange coupling 
J, and in La2Cu04 or Sr2Cu02Cl2 (2122) this is about 
J ~ 130 meV. Although thermal neutron scattering is 
one of the most powerful tools in probing excitations in 
magnetic systems, it is not normally suitable for study- 
ing the physics at such a high energy scale. As a re- 
sult, neutron scattering studies of the long- wavelength 
spin dynamics of 5 = 1/2 2DQHA's such as La2Cu04 
are very limited . The recent discovery of the sys- 
tem Sr2Cu304Cl2 (2342), which contains a square lattice 
Cu// subsystem that has an order of magnitude smaller 
J, has opened the door for a quantitative study of the 
spin dynamics of the 5=1/2 2DQHA @. In this 
Letter, we report systematic measurements of ujq as a 
function of T/J in the model 5=1/2 2DQHA's, 2122 
and 2342. We show that dynamic scaling is obeyed over 
the entire temperature range probed in our experiment, 
0.2 < T/J < 0.7, which corresponds to the correlation 
length range 100 > £/a > 2. 

The neutron scattering experiments were carried out 
at the Center for Neutron Research at the National In- 
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stitute of Standards and Technology. Large single crys- 
tal samples of both 2122 and 2342 were grown by the 
flux method; typical sample dimensions were 20 x 20 x 5 
mm 3 , and the sample mosaicity was ~ 0.2° half width at 
half maximum (HWHM). Various experimental configu- 
rations at the thermal beamlinc BT9 were employed in 
the measurements. Progressively coarser resolution was 
used to accommodate the wider energy width as the tem- 
perature was increased and, correspondingly, £ decreased. 
Pyrolytic graphite filters were used to reduce the effects 
of A/2 neutrons. The magnetic structure factor at (71", n) 
as a function of energy transfer, cj, was measured at var- 
ious temperatures. The lower bound of the temperature 
range was governed by the Neel ordering temperature 
(TV « 256 K for 2122 and T N « 40 K for 2342), be- 
low which the spin dynamics is no longer relaxational at 
(tt, 7r) pUTl]]. Typical scans are shown in Fig. |l| for both 
2122 and 2342. The energy resolution is also plotted in 
each panel MS M. horizontal bar. 

To determine ujq these scans were fitted to the following 
form for the dynamic structure factor convoluted with 
the instrumental resolution function: 
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Here q is defined as the deviation in wave vector from 
the 2D antiferromagnetic position (7r,7r), and c is the 
spin wave velocity. We have used T q = wq[1 + ^{qC) 2 ] 1 ^ 2 
as proposed by Tyc, Halperin, and Chakravarty (THC) 
Our fits turned out to be insensitive to the choice 
of the phenomcnological parameter fx, and, accordingly 
/1 was set equal to zero. In Eq. (Q), we define the char- 
acteristic energy scale as the HWHM of the quasielastic 
peak, that is, luq = r q= o. It should be noted that Eq. 
(Q) is a simplified version of the dynamic structure factor 
obtained by THC from their molecular dynamics simu- 
lation of a classical lattice rotor model |l2] . Specifically, 
we ignore the logarithmic corrections (log[l + (q£) 2 ]) in 
the dynamic structure factor of THC, since we only focus 
on the small q regime. 

We plot the results for uj^j J versus the inverse corre- 
lation length £ _1 in Fig. |2[ The values of the correlation 
length in these systems have been determined previously 
with high accuracy ||[ll|]. Therefore, ujq and the overall 
amplitude 5(0) are the only adjustable parameters in the 
fits to Eq. (Q). In order to compare the two different sys- 
tems, as well as with the quantum Monte Carlo results, 
we scale loq by J, and £ by the lattice constant a. Each 
data set obtained with a different experimental setup is 
plotted in a different symbol. Filled and open symbols 
are used for 2122 and 2342, respectively. In their study 
of the spin dynamics of La2Cu04, Hayden and coworkers 
obtained u; at T = 320 K M, which is also plotted in 



Fig. 0. It is quite remarkable that the experimental re- 
sults from all three systems fall on a single straight line 
without any adjustable parameter. This clearly shows 
the dynamic scaling behavior of loq. If we fit the data 
to the scaling form loq ~ £~ z , we obtain z = 1.0(1), 
in excellent agreement with the theoretical value for the 
2DQHA. The solid line is a fit of the experimental data 
to the scaling form with fixed z = 1: cua/J = l.l(l)a£ _1 . 

In order to compare our experimental results with 
the results of numerical simulations, Makivic and Jar- 
rell's quantum Monte Carlo results jl3| for the 5=1/2 
2DQHA are plotted in Fig. ||. Wysin and Bishop's Monte 
Carlo molecular dynamics calculation results for the clas- 
sical 2D Heisenberg antiferromagnet are also plotted jl4} . 
Not surprisingly, the quantum Monte Carlo results are 
in an excellent agreement with our neutron scattering 
results. However, it is interesting to note that the val- 
ues from the classical Monte Carlo calculation also agree 
with our neutron scattering results within experimental 
error bars. Since Wysin and Bishop have calculated £ 
as well as loq, there is no adjustable parameter except 
the J — > J5(5 + 1) scaling. The good agreement be- 
tween quantum and classical dynamics presumably re- 
flects the fact that the zero temperature spin dynamics 
of the 5 = 1/2 2DQHA are well described by the classical 
spin- wave picture with a uniform frequency renormaliza- 
tion of only 17%. 

Two important comments regarding the data analysis 
are in order. First, due to the nonzero resolution width in 
both energy and momentum transfer, scattered neutrons 
with q 7^ are invariably included in our data shown in 
Fig. 1. However, from numerical simulations, we have 
verified that the q / dynamical contribution to the 
total quasielastic scattering intensity for 2342 is insignif- 
icant and yields a correction to the energy width that 
is much smaller than the experimental error bars. On 
the other hand, 2122 has an order of magnitude larger 
spin wave velocity (820 meVA) than that of 2342 (95 
meVA), and accordingly 2122 has a very steep disper- 
sion relation. For 2122, the nonzero q dynamical contri- 
bution thus gives an apparent peak width in energy that 
is slightly larger than the intrinsic value. Although this is 
already taken into account in our fitting processes, rela- 
tively large error bars are generated for 2122 as a result. 
We have employed a number of different experimental 
configurations to ensure that the extracted ujq is intrin- 
sic in both 2122 and 2342. Second, we have defined the 
characteristic energy scale as the HWHM of the dynamic 
structure factor at q = 0. In Monte Carlo studies of Ref. 
Jl3[ | and Ref. [H, a different definition of loq was used, 
and the values given in these studies are converted to 
fit our definition by multiplying by a constant. We have 
obtained this constant multiplication factor of 0.682 by 
fitting the raw data in Ref. @ to Eq. (Q). 

The consistency among the experimental results and 
the numerical simulation results in absolute units gives 
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strong credence to the combined data as a testing ground, 
against which various analytic theories can be examined. 
Since the pioneering work by CHN, there have been a 
number of theoretical studies to determine the charac- 
teristic energy scale of the 2DQHA. They all seem to 
indicate z = 1 dynamic scaling behavior of luq, but dif- 
fer in the temperature dependence of R^, which is de- 
fined as a dimcnsionless ratio = wo£/co- Specifically, 
Ru> ~ T 1 ' 2 is predicted by CHN, while Grempel |J| has 
carried out a conventional mode-mode coupling calcula- 
tion and obtained the explicit prefactor of the T 1 / 2 de- 
pendence: 
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This result is close to the expression obtained by THC 
in their numerical simulation: R w — 0.85y // T/2'Kp s w 
0.8y/T/J. However, Auerbach and Arovas, in their 
Schwinger boson mean field theory |TJ], have obtained 
a T-independent R w . Quantum Monte Carlo calcu- 
lations by Makivic and Jarrell |Tj| as well as Nagao 
and Igarashi's self-consistent theory |0 also show a T- 
independent although these results are limited to a 
somewhat narrow temperature range. 

In Fig. |5[ we plot the dimensionless ratio i? w obtained 
from our measurements as a function of reduced tempera- 
ture. Note that cq for the S =1/2 2DQHA is well known 
from both Monte Carlo |Q and series expansion stud- 
ies ||l8| : co = 1. 657 Ja. Therefore, Fig. ^ also does not 
contain any adjustable parameter. Results from Monte 
Carlo calculations are also plotted in the figure, where 
the same symbols as in Fig. |^ are used. The analytic 
theory by Grempel, Eq. (||), is plotted as a solid line, 
while a fit to a presumed constant R u is plotted as a 
dashed line. 

We should emphasize that the dominant T dependence 
of wo is ~ £ which has a strong exponential T depen- 
dence, so that Fig. |^ only shows the weak T dependence 
of the prefactor of the leading £ -1 behavior. The salient 
feature to recognize in this figure is that the coupled 
mode theory of Grempel is in a reasonable agreement 
with our experimental results in absolute units. Consid- 
ering that the mode-mode coupling approximation used 
by Grempel ignores vertex renormalizations, the agree- 
ment in this intermediate temperature range appears to 
be quite good. The large error bars on our 2122 data, 
however, prevent us from concluding that R u indeed has 
the predicted y/T dependence. A temperature indepen- 
dent Ruj describes our data equally well, if not better. 
Another analytic theory that predicts a T-independent 
R u is quantum critical theory |19|. However, quantita- 
tive comparison with our data is difficult, since the tem- 
perature dependence of loq in the quantum critical regime 
is not known theoretically. 

It turns out that the task of experimentally determin- 



ing the functional form of R u is very difficult. Since vT 
is a weak function of T, one has to probe at very low 
temperatures to distinguish between \[T and constant 
behavior. Although Sr 2 Cu0 2 Cl 2 has the smallest Tj^/J 
among model S = 1/2 2DQHAs known to this date, its 
large nearest neighbor coupling makes a high-resolution 
neutron scattering study very difficult. To elucidate fur- 
ther the temperature dependence of the correction term 
R u , a quantum Monte Carlo study of the spin dynam- 
ics over a wide temperature range is highly desirable. It 
should be noted that in the bicritical S = 5/2 2DQHA 
system Rb2MnF4, R w appears to depend strongly on T; 
indeed in that case the observed temperature dependence 
is much more rapid than T 1 / 2 [ ppf . 

In summary, a neutron scattering study of the long- 
wavelength dynamic spin correlations in the model S — 
1/2 2DQHA's Sr 2 Cu0 2 Cl 2 and Sr 2 Cu 3 4 Cl 2 has been 
presented. We have measured the characteristic energy 
scale as a function of temperature, and have shown that 
dynamic scaling is valid over a wide temperature range 
(0.2 < T/J < 0.7) with the dynamic critical exponent 
2 = 1. This temperature range corresponds to the mag- 
netic correlation length of 100 > £/a > 2, for which 
dynamic scaling is expected to hold according to theo- 
retical predictions. Our attempt to determine the weak 
temperature dependence of any possible corrections to 
this z = 1 scaling is, however, inconclusive, and remains 
to be addressed in future studies. 
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FIG. 1. Representative a;-scans at the position, 
which corresponds to (1/2 1/2 0) for both Sr2Cu304Ci2 and 
Sr2Cu02Cl2, following the notation of Refs. [3,11]. The back- 
ground from incoherent scattering, measured far away from 
the (n,n) position, has been subtracted from the raw data. 
A representative scan of the incoherent scattering, multiplied 
by 1/5, is shown as a dashed line in the top right panel. A 
fixed final neutron energy of 30.5 meV and collimations of 
40'-48'-Sample-44'-40' have been used. The solid lines are 
the results of least square fits to Eq. ([!]) convoluted with the 
instrumental resolution function. 
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FIG. 2. Logarithmic plot of the reduced characteristic en- 
ergy scale loqI J versus the inverse of reduced spin correlation 
length a/£. The open symbols are the data for Sr2Cu304Ci2 
plotted with J — 10.5 meV, and the filled symbols are the 
data for Sr 2 Cu0 2 Cl2 plotted with J = 125 meV. Note that 
we use different symbols for each experimental configuration. 
Here HKO denotes the scattering plane in reciprocal space. 
We also plot the result of Ref. [9] for La 2 Cu0 4 at T = 320 K. 
The solid line is the dynamic scaling prediction with z = 1, 
that is, u>o ~ The results of the quantum Monte Carlo 

simulations of the S — 1/2 2DQHA by Makivic and Jarrell 
[13] are plotted as *, while Wysin and Bishop's classical Monte 
Caro result [14] are plotted as x with the energy scaled by 

js(s+i). 
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FIG. 3. The dimensionless ratio = luq^/cq as a function 
of reduced temperature. The same symbols as in Fig. ^ are 
used in the plot. The solid line is the result from Grempel's 
mode- mode coupling theory, Eq. (^). The dashed line is a fit 
to a constant, — 0.66(7). 
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